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Abstract. We develop two iteration scliemes for construction of localized stationary solutions 
(bumps) of a one-population Wilson-Cowan model with a smoothed Heaviside firing rate func- 
tion. The first scheme is based on the fixed point formulation of the stationary Wilson-Cowan 
model. The second one is formulated in terms of the excitation width of a bump. Using the 
theory of monotone operators in ordered Banach spaces we justify convergence of both iteration 
schemes. 



Neural field models have been the subject of mathematical attention since the publications 
mill SHI. These models typically take the form of integro-differential equations. We consider 
a one-population neural field model of the Wilson-Cowan type HI ID [3l IH |5l 

-|-oo 



Here u{x, t) represents the activity of population, / the firing-rate function, co the connectivity 
function, and h the firing threshold. For review on the model ( |1.1| ) see [5]. Existence and stabil- 
ity of spatially localized solutions and traveling waves are commonly studied for the case when 
the firing rate function is given by the unit step function JHOm. However, the results for the 
case when the firing rate function is smooth are few and far between fTllSl fTOlfTTI . 

In the mathematical neuroscience community time-independent spatially localized solutions of 
( |1.1[ ) are referred to as bumps. The motivation for studying bumps stems from the fact that they 
are believed to be linked to the mechanisms of a short memory [12|. In the case when / is 
given as a unit step function, one can find analytical expressions for the bump solutions H. 
In principle, bumps solutions can also be constructed when the firing rate function is smooth 
provided the Fourier-transform of the connectivity function is a real rational function. In that 
case the model can be converted to a higher order nonlinear differential equation which can 
be represented as a Hamiltonian system. The bumps are represented then by homoclinic orbits 
within the framework of these systems. See for example |[T3l[T4llT5l . 

Kishimoto and Amari Q have proved the existence of bump solutions of ( |1.1[ ) when / is a 
smooth function of a special type (smoothed Heaviside function), using the Schauder fixed point 
theorem. The Schauder fixed point theorem, however, does not give a method for construction of 
the bumps. Pinto and Ermentrout in [9J constructed bumps using singular perturbation analysis. 
However, this method is quite involved, and is restricted to the lateral-inhibitory connectivity 
(i.e., io is assumed to be continuous, integrable and even, with a;(0) > and exactly one posi- 
tive zero). Coombes and Schmidt in |8| developed an iteration scheme for constructing bumps 
of the model ( |1.1| ) with a smoothed Heaviside function. They, however, did not give a mathe- 
matical verification of their approach. Apart from the work of Coombes and Schmidt [81, the 
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authors of the present paper do not know about other attempts to develop iterative algorithms 
for the construction of bumps. Thus there is a need for a more rigorous analysis of iteration 
schemes for bumps. This serves as a motivation for the present work. 

We present two different iteration schemes for constructing bumps. The first one is based on the 
fixed point problem introduced in |7]. The second scheme, which is modification of the proce- 
dure introduced in HI, is an iteration scheme for the excitation width of the bumps. We prove 
that both schemes converge using the theory of monotone operators in ordered Banach spaces. 

The present paper is organized in the following way: In Section |2] the properties of the one- 
population Wilson-Cowan model are reviewed with emphasis on the results of Kishimoto and 
Amari fT]. In Section |3] some necessary mathematical preliminaries are introduced. Section [4] 
is devoted to the study of a direct iteration scheme based on the fixed point problem proposed 
by of Kishimoto and Amari Q. In Section |4~T] we illustrate the results with a numerical exam- 



ple. In Section [5] we introduce a fixed problem based on the specific representation of the firing 
rate function studied in [8l. The fixed problem is formulated for the crossing between bumps 
and a shifted parameterized threshold value h + t, t > 0. The bump solution can be restored 
from these crossings. We prove that there is a fixed point which can be obtained by iterations. 
We provide an numerical example in Section 5.1 In Section [6] we summarize our findings and 
describe open problems. 



2. Model 

Let / : M — [0, 1] be an arbitrary non-decreasing function. We assume that the connectivity 
function uj satisfies the following conditions: 

(i) UJ is symmetric, i.e. uj{—x) = uj{x), 

(ii) f^\u{x)\dx < CO, i.e., UJ e L^{R), 

(iii) UJ is continuous and bounded, i .e., uj G i?C(M), 

(iv) UJ is differentiable a.e. with bounded derivatives, i.e., uj G 1/F^'°°(M). 

The examples of such a function are 

(2.1) uj{x) = A'e"''^" - Me-""^", < M < K, < m < k, 
and 

(2.2) a;(2;) = e"^l^l(6sin|x| +cosx), 6 > 0. 



In Fig{TJa) we illustrate the function given in (2J_i with parameters K = 1.5, k = 2, and 



M = m = 1. In FigjTJb) we illustrate the function in (2/2) with b = 0.3. The function ( 2.1 
models the lateral-inhibition coupling and is often called as a Mexican-hat function, e.g., see 
||4l[l4l|5l. The model with periodically modulated spatial connectivity given by ( |2.2[ ) was con 
sidered in 



Stationary solutions of ( |1.1| ) are given as solutions to the integral equation 

+ 00 

(2.3) u{x) = / w{y - x)f{u{y) - h)dy. 



We note the following properties of (2.3 1 



A solution u is translation invariant. That is, if u{x) is a solution to (23 1, so is u{x — c) 
for arbitrary constant c S M. 



ITERATIVE SCHEMES FOR BUMP SOLUTIONS IN A NEURAL FIELD MODEL 



3 



0.5 










0.4 






0.3 






0.2 






0.1 

















-0.1 






-0.2 







1 

0.8 
0.6 
0.4 
0.2 



-0.2 



-0.4 





-20 



-10 



10 



20 



(a) 



(b) 



Figure 1. Examples of the connectivity function uj: (a) The Mexican-hat 
function \2A\ , and (b) the function \2.2\ , with the parameters given in the text. 



A symmetric solution to (2.3 1 can be expressed as 



(2.4) 



u{x) = J r{x,y)f{u{y) - h)dy 





where 



r{x, y) = w{y - x) + w{y + x). 
Let the function / be given as the unit step function 



(2.5) 



f = e, 9{u) 



0, M < 

1, u>0. 



Amari was the first who observed that in this case, the spatially localized solutions to ( 2.3 1 
can be explicitly constructed. Following [4 | we introduce the following definitions: 

Definition 2.1. The set R[u] = {x\ u{x) > h} is called the excited region ofu{x), lH. 

Definition 2.2. An equilibrium solution u{x) of ( |1.1[ ) with f = 6 is called a bump with the 
width a, if the excited region of u is an interval of the length a, i.e., i?[n] = (01,02), where 
a = 02 — ai. 

Then a bump solution with the width a is given as 

u{x) = J uj{y — x)dy, a = 02 — ai. 

ai 



Due to translation invariance of ( |2.3| ) we without loss of generality consider bumps defined on a 
symmetric interval, i.e., 

a/2 



u{x) = J u}{y — x)dy, 02 = — ai = a/2. 

-a/2 
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It is easy to see that u{x) in this form is a symmetric function. Indeed, letting z = —y we have 

-a/2 a/2 
a/2 -a/2 



Thus, using ( |2.4| ) a bump solution can be written as 

a/2 

(2.6) u{x) = j r{x,y)dy. 



We define a new function <I> 

ry 

^{x,y) = / r(x, X, y G M, y > 
Jo 

with 

(2.7) — (x,y) =tj(y + x) -tj(y-x), — (x, y) = r(x, y). 
We conveniently express bumps by means of the function 

Theorem 2.3. Let h > be fixed. The model ( |1.1[ ) vviY/z the firing-rate fimction f = 9 possesses 
a bump solution if and only if there exist a width, a, such that 

(2.8) <^{a/2,a/2)= [ w{y)dy = h 

Jo 

and 

(i) $(x,a/2) < /i, Vx > a/2, 

(ii) $(x,a/2) > /i, Vx e [0,a/2). 

The bump solution is given then as u{x) = <I>(x, a/2). 

The stability of bumps has been studied using the Amari approach [4] and the Evans function 
technique, [5|. Here we present the result based on [4|: 

Theorem 2.4. Let h > be fixed, f = 9, and there exist a bump with the width a. The bump is 
linearly stable ifuj{a) < and unstable ifuj{a) > 0. 

The firing-rate function we treat here is of the following type, Q 

To, n < 

(2.9) /(n) = <^ ^(u), < n < r , 

[ 1, n > r 

where r > 0, i;Ms an arbitrary continuous, monotonically increasing, and normalized function 
such that 

m = 0, <p{t) = 1. 

The example of such a function is 

' 0, n < 



(2.10) /(u) = S(n/r,p), 



nP + (1 - u)P 
1, u> 1 



, <u<l , p> 0, 



where S(-,p) € C^^^M) and [p] denotes the integer part of p. 



We need the following definition: 
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Definition 2.5. R* [u] = {x\u{x) > h + t} is called a maximally excited region, and R [u] = 

{x\h < u{x) < h + t} is an incompletely excited region, ff]. 

Definition 2.6. An equilibrium solution u{x) of with f given by (2.9 1 is called a bump if 
R* [u] is the interval surrounded by an incompletely excited region R^ [u] , i.e., R[u] = R* [u] U 
R^ [u] being another interval, |7|. 

Thus, by Definition 2.6 the function u{x) displayed graphically in Fig|2]can be a bump to 
with / given in (2^1 when r = ti, whereas for r = T2 it can not be a bump. 




Figure 2. The graph of a function u{x) which satisfies Definition 2.6 when 
T = Ti and which does not satisfy it when r 



T2- 



Let / be given as ( |2.9| ), fo{u) = 9{u), and fr{u) = 0{u — r). To distinguish between bump 
solutions to ( |1.1| ) with different firing rate functions, we use the following terminology: the 
neural field with the firing rate functions /o, /r, and / is called a /o-field, /T--field, and /-field, 
respectively. We observe that /T--field is equivalent to the /o-field with the new threshold value 

/i + r, and 

fr{u) < f{u) < /o(n). 

The original idea of Kishimoto and Amari [7 | is to use bump solutions of the /q- and /T--fields 
to prove the existence (and stability) of bumps in the /-field. If lo has a Mexican-hat shape (e.g., 
see Figjlja)) then the /o-field (/r-field) possesses two symmetric bumps for moderate values 
of h, one stable and one unstable bump. In [7| it was shown, using the Schauder fixed point 
theorem, that there exists a bump solution of /-field if both /q- and /T--fields possess linearly 
stable bumps and uj has a Mexican-hat shape (i.e., the connectivity function can have the shape 
as in Fig[T];a) but not as in Fig|TJb)). Moreover, if is a differentiable function it was shown 
that the /-field bump is stable. 

Notice that the differentiability of (p can be replaced by a weaker assumption, namely differ- 
entiability almost everywhere, i.e., (p G 1^^'^([0, r]). Then, the firing-rate function ( |2.9[ ) can be 
represented as in |[8l, i.e., 



(2.11) 



with 6 given by ( |2.5| ), suppjp} = [0, r], and p is positive and normalized J p{x)dx = 1. 

— oo 

In this paper we prove existence of bumps in the /-field, and introduce two iteration meth- 
ods for their construction. We improve the existence result obtained in Q by relaxing on the 
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assumption that a; has a Mexican-hat shape. We also do not require the bumps of the /g- and 
/T--fieId be stable as it is assumed in [7|. So far there have been two methods used to construct 
bumps in /-field: One is based on the singular perturbation analysis, [9|. This method is quite 
involved and, moreover, it restricts the choice of a; to functions of a Mexican-hat shape. The 



other method is to convert ( |2.3| ) to a higher order nonlinear differential equations which can be 
represented as a Hamiltonian system. The bumps then are given by homoclinic orbits within 
the framework of these systems, see |[T3l[T4l[T5l . This method requires the Fourier transform 
of to be a real rational function. Thus, it can not be applied in some cases, as for example in 



the case of ( |2.1| ). We do not requite uj to have either Mexican-hat shape or real rational Fourier 
transform to be able to apply our iteration schemes. 
We use the following assumptions: 

Assumption 1. There exists a bump with the width 2Ao of the fo-field model, and a bump with 
the width 2 At- to the fr-field model. Moreover, the widths are such that A,- < Aq. 

To illustrate this assumption let us assume that there is a bump solution of the /o-field model, 



i.e., $(Ao,Ao) = h, see Theorem 2.3 Then, by the inverse function theorem there exists 
a value r > such that <&(At-, A,-) = /i + r for some A,- < Aq, if c<j(2A) < in some 
vicinity of Aq . In this case both bumps are stable by Theorem |2.4[ However, Assumption [T] 
can be satisfied even when the situation described above does not take place, i.e., the condition 
a;(2A) < is not fulfilled for all AG [A^, Aq], see for example Figjs] 

Under Assumption[T]bumps for the /o-field model and the /r-field model are, in accordance 
with Theorem |2. 3 1 given as 

uo{x) = Ao) Ur{x) = ^{x, A^). 
Assumption 2. The function r{x, y) is non-negative for all x,y £ [A,-, Aq]. 

We get the following relationship between Ur and uq: 
Lemma 2.7. Under As sumption^we have Ur < uq on [A,-, Aq]. 
Proof. We get 

/•Ao /-Ar /"Ao 

UQ{x)-Ur{x)= I r{x,y)dy- / r{x,y)dy = / r{x,y)dy >0. 



JAr 

□ 

In this paper we will only consider bump solutions of the /-field such that 
(2.12) u{x) > h + T, \/x £ R[ur -t], u{x) < h\/x ^ R[uo\. 

3. Mathematical Preliminaries 

Let K be. a. cone in a real Banach space B and < be a partial ordering defined by K. Let 
w,v G ;B be such that w < v. Then a set of all g £ B such that w < g < v, defines an ordered 
interval which we denote [tt', w]] • 

The theoretical foundation of the iteration schemes presented in Section]?] and Section ]5] are 
based on the following general results: 

Theorem 3.1. Let wo,vo G B, wq < vq and A : [twc^ol ~^ B be an increasing operator 
(Aw < Av provided w < vfor any w,v £ B) such that 

wq < Awo, Avq < Vq. 

Suppose that one of the following two conditions is satisfied: 
(HI) K is normal and A is condensing; 

(H2) K is regular and A is semicontinuous, i.e., — )• x strongly implies Axn — )• Ax weakly. 
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Then A has a maximal fixed point x* and a minimal fixed point x=k in \wq, vq^ ; moreover 



X* = lim Vn, X* = lim Wn, 



where Vn = Avn-i and Wn = Awn-i, n = 1,2, 3..., and 

Wo < Wi < ... <Wn< ••• < Vn < ••• < < Vq. 

See wn . 



From Theorem 3.1 we get the following corollary. 



Corollary 3.2. If under the conditions of Theorem 3.1 then X is the unique fixed 

point of the operator A in \wq, vq^. 

Theorem 3.3. The cone K = {u £ B\u{x) > 0} is normal but not regular in B = C{D), and 
regular in B = Lp{D), \ < p < oo, where D is a bounded set and D is a closed bounded set, 
see wn . 

Theorem 3.4. The Hammerstein operator 

{Af){x)= t k{x,y)ij{y,f{y))dy 



is continuous and compact in C{[a, h]) ifk{x, y) and'ip{x, y) are continuous functions on [a, b] x 
[a,b]. 

Proof. The operator A can be represented as the superposition, A = LN, where L is the hnear 
operator 

{Lg){x) = / k{x,y)g{y)dy, 

J a 

and N is the Nemytskii operator 

{Nf){x) = i^{x,f{x)). 

The linear operator L : C{[a,b]) — M is continuous and compact if k{x,y) is continuous 
|fT9i l. Obviously, the Nemytskii operator N : C{[a, b]) — C{[a, b]) is continuous and bounded 
if ip{x,y) is continuous. Thus, the Hammerstein operator A is completely continuous as the 
superposition of the continuous and bounded operator N, and completely continuous operator 
L. □ 



4. Iteration Scheme I: Direct Iteration. 

In this section we consider the direct iteration scheme for construction of bumps. This scheme 
is based on fT|. We start out by observing that a bump solution of an /-field satisfying ( |2.12 1 
can be written as 

(4.1) u{x) = Ur{x)+l r{x,y)f{u{y) - h)dy, 

for all X e [A^, Ao]. 

First, we prove that there exists a solution u* {x) to ( |4.1| ) and it can be iteratively constructed. 
Next, we introduce assumptions under which u* {x) is appeared to be a restriction of a bump 
solution to ( |1.1[ ) on [A,-, Aq]. Finally, in Section we illustrate our results numerically and draw 
some conclusions based on the numerical observation. 



Let S be a real Banach space with partial ordering > defined by the cone K = {u £ 
B\u{x) > 0}. We have the following theorem. 
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Theorem 4.1. LetB be either L2{[Ar, Aq]) or C{[Ar, Aq]). Let oj satisfy As sumption^and^ 
the operator Tf : , uq} C B ^ B be defined as 

/■Ao 

(4.2) {Tfu){x) =Ur{x)+ I r{x,y)f{u{y) - h)dy. 

JAr 

Then Tf has a fixed point in \ut-, uq^. Moreover, the sequences {Tjur} and {TJuq} converge 
to the minimal and maximal fixed point of the operator Tf, respectively. 



Proof. We base our proof on Theorem |3.1| The cone K is normal provided B = C{[Ar, Aq]) 
and is regular provided B = L2([A,-, Aq]), see Theorem 3.3 By Assumptions [T] and |2] there 
exist Ut and uo such that < Ur{x) < uo{x) for all x G [A,-, Aq] . Thus, {ur, uq} is the ordered 
interval defined on B. 

We describe the properties of Tf which hold true in both spaces: B = C([At-, Aq]) and B = 
L2([A,-, Aq]). First of all, Tf is positive and monotone due to Assumption |2] and monotonicity 
of /, i.e., 

Ul(x) < U2{x) =^ {TfUl){x) < (TfU2)ix). 

Moreover, Tf is continuous due to continuity of / and boundedness of r{x, y). Defining a non- 
linear operator Tg associated with the non-negative function g{x) by 

/■Ao 

{Tgu){x) = Ur{x) + I r{x,y)g{u{y) -h)dy 

JAr 

we get 

T/gno = no, Tf^Ur = Ur- 
Thus, from Assumption [2] we can easily deduce that 

g{x) < m{x) =^ {Tgu){x) < (Tmu){x), 

and, therefore, 

(T/,n)(x) < {Tfu){x) < {Tf,u){x). 

We obtain 

TfUr > Tf^Ur = Ur, TfUQ < Tf^UQ = Uq. 

From Theorem 3.1 we conclude that Tf : lur, uq} C L2([At-, Aq]) — )■ L2{[Ar, Aq]) has a 
fixed point in lur, uo|| which can be found by iterations. However, for the case S = C([At-, Aq]) 
it remains to show that Tf is condensing. Applying Theorem 3.4 to the Hammerstein operator 

on the right hand side of ( |4.2[ ), i.e., to the operator u i— > r{x, y)f{u{y) — h)dy, we find that 
Tf : lur,uo} — )• C([A,-, Ao]) is compact and, thus, condensing. This observation completes 
the proof. □ 

Next we show that the fixed point u* of the operator Tf refened to in the theorem above can 
be extended to the solution u of (23 1 over M in such a way that u{x) > /i+r for x G [0, A,-] and 
u{x) < /i for X G [Aq, oo). To do so we introduce additional assumptions on the connectivity 
function uj. 

Assumption 3. uq is a decreasing function on the interval [Ar, Aq] which is equivalent to 

— (x,Ao) < 0, Vx G [A^,Ao], 
and Ur is a decreasing function on [Ar, Aq] which is equivalent to 

— {x,Ar) < 0, Vx G [A^, Ao]. 
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From this assumption the transversahty of the intersections uo{x) with h, and Ut-{x) with 
h + T follows. Thus, the assumption always can be satisfied if, for example, we choose a small 
r provided | Aq — A,- 1 is sufficiently small. 

Assumption 4. 

/ I T^{x,y) I /(uo(y) - h)dy < - — (x,A^), Vx G [A^, Aq]. 

Assumption[4]is technical and is used to prove that u* (x) is a decreasing function on [A,- , Aq] . 
Noticing that 

dr _ 

dx dxdy ' 

non-negativity of dr/dx for x,y G [A,-, Aq] and Assumption 3 imply that Assumption 4 is 
satisfied. Indeed, the following chain of inequalities is valid for all x E [A,-, Aq] 

However, the non-negativity of dr / dx is rather rigid condition. 



Lemma 4.2. Let Assumption 3 and 4 be fulfilled. Then, the fixed point u*{x) of the operator Tf 
is differentiable and decreasing on the interval [A,- , Aq] . 



K(x))'=<(x)+I, /= / {x,y)f{u*{y)-h)dy. 



Proof. We get 

f Ao 

JAr dx 

In order to prove that (u*(x))' < we need to show that / < — n'^(x) where 

<(x) = ^(x,A,) < 
ox 

by Assumption [3] 

We get the following chain of inequalities for |I| 

\I\ < J^' \ ^{x,y) \ f{u*{y) - h)dy < 
/•Ao 

< I Q^(^^y) I /("o(y) - h)dy. 

Thus, by Assumption |4] we have |/| < — n'^(x) and therefore / < —u'^{x). □ 

Finally, we introduce the assumption which by Definition 2.6 allows us to view the extended 
solution u of u* as a bump: 

Assumption 5. The function $ is such that 

(i) ^{x,y) < /i, Vx > Ao, y £ [A^, Aq], 

(ii) «>(x, y) > /i + r, Vx e [0, A^], y € [A^, Aq]. 

Theorem 4.3. Let u* G C^([At-, Aq]) define the fixed point of the operator Tj referred to in 
Theorem \4.1\ Under Assumptions^-^the function u* can be extended to a bump solution u{x) 
of ( |2.4[ ) defined on M in such a way that n(x) > h + t for all x G [0, A,-) and u{x) < hfor all 
X G (Aq, 00). 
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Proof. From Theorem 4.1 and Lemma 4.2 it follows that there exist unique 5r, 5o ■ A,- < 5r < 
6o < Aq such that 

u*{6r) = h + T, u*{6o) = h. 
Let us introduce the function F defined by 

1, 0<y<6r 
Fiy) = { f{u*{y)-h), 6r<y<6o 
0, y> 6o. 



Then, according to Lemma 4.2 F{y) is monotonically decreasing function on [^n^o] with 



F{6o) = and F{6r) = 1. From (|23]l we get 

rSo rSo rSo 

u{x) = I r{x,y)F{y)dy = — / r{x,y) / F'{z)dz. 
Jo Jo Jy 

By changing the order of integration, we have 

/•So j-z i-So j-z 

u{x) = - / r{x,y)dyF'{z)dz = / r{x,y)dyd{l - F{z)), 

JSr Jo JSr Jo 

or 

u{x) = / / r{x,y)dydS, = / ^{x, z{(^))d^, where ^ = 1 — 
Jo Jo Jo 
It remains to show that <I>(x, z) < hfor x > Aq, z E [6r, Sq], and ^{x, z) > /i + r for < x < 
A^, z & [Sr, 5o]- Assumption [5] guarantees that these inequalities are fulfilled even on a larger 
interval. Moreover, u{—x) = u{x) due to symmetry of uj. Thus, the proof is completed. □ 



The proof of Theorem |4.1| is a modification of the theorem used in [7]. The modification is 
caused by the fact that our assumptions on the connectivity functions uj are different from ones 
used in fl\. 

4.1. Numerical example. In this section we exploit examples of a; be given by ( |2.1[ ). Thus, 
the equation uj{x) = has one positive solution Xmax = \/l/(A: — m) hi{K/M). Furthermore, 
oj{x) is positive for all \x\ < Xmax and is negative for all |x| > Xmax- This defines the behavior 
of the antiderivative to uj, i.e., 

X 

(4.3) Wix) = J uj{y)dy. 



Then, for any /i, r satisfying the inequality 

lim W{x) < h < h + T < W{Xmax) 

there exist the widths A,-, Aq referred to in Assumption 1 such that [A^, Aq] C {l/2xmax,oo.) 
Moreover, ^0(2;) = ^{x, Aq) and Ur = ^{x, A,-) are bumps of the /q- and the /T--field model. 

We let K = 1.5, k = 2, M = m = l(see Fig{T]), and chose h = 0.1 and r = 0.05. Given 



these parameters and ( |2.8[ ) we found Aq = 0.6633 and two possible values of A,- : 0.1769 and 
0.5012 (see Figjsja)). We fix A^ = 0.1769, and denote A^* = 0.5012. We notice here that 



a;(2Ao) and uj{2Af^) are negative while a;(2AT-) is positive, see Fig. [Sj^b). By Theorem 2.4 
we conclude that uq{x) = <I>(x, Aq) and n^* = ^{x, A**) are linearly stable solution to /q- and 
/T--field models, while Ur{x) = ^{x, A,-) is a linearly unstable solution to /T--field model. This 
explains upper index 'st' in our notation. 



Assumptions 2-5 has been verified numerically. Thus, we apply Theorem 4. 1 We chose / 
to be given as in in ( |2.10| ) with p = 3. In FigQa) we have plotted the fixed point u* {x) of the 
operator Tf obtained by iterations from the restriction of Ur{x) and uo{x) on [A,-, Aq] as the 
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Figure 3 . (a) The antiderivative of 00 given by \A.2>) . The blue square indicates 
the point (2At-, /i + r), the blue circle corresponds to (2A^*,/i + r), and the red 
circle corresponds to (Aq*, h). The horizontal solid and dashed line are h and 
h + T, respectively, (b) The connectivity function and the points (o, cj(a)) 
with a = At, A^* and Aq, denoted by the blue square, blue circle, and red 
circle, respectively. 



initial values. From Corollary 3.2 we conclude that u* is, then, a unique solution of the fixed 
point problem for ( |4.2[ ). We have plotted uq, Ur, and u^* on the same figure to illustrate the 
inequality 

(4.4) <* <u* <uo. 

Thus, u* is located in between of two stable bumps of the /q— and the fr field model on 
[A,-, Aq]. Based on [4| we claim that u* is a restriction of the stable bump solution to /—field 
equation. 

FigQb) illustrates the dynamics of the iteration process. There we have plotted the numerical 
errors calculated as 

(4.5) ein) = max\{Tfuo){x)-{T^Ur){x)\, n = l,2,...,N, 

where Tj? is equivalent to the identity operator, n corresponds to the iteration number, and A'^ 
denotes the total number of iterations. We observe that in our calculations e{n) < 10^^ for 
n > 16. 

In Figj5]we have plotted the graph of the bump solution to the /-field model obtained as the 



extension of u* from [A,-, Aq] to M, see Theorem 4.3 



Assumptions 1-5 are satisfied for some set of parameters when co is given as in ( |2.2[ ). We, 
however, do not provide this example here. 



5. Iteration Scheme II: Bumps Width Iteration 

In m an iteration procedure for construction of bump solutions to the /-field model has been 
worked out. However, a mathematical verification of the procedure has not been given. In the 
present section we introduce an iteration scheme which is based on the idea presented in [8J 
and give a mathematical verification of this approach, we illustrate the results with a numerical 



example using the same parameters as in Section 4. 1 
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Figure 4. (a) A fixed point of the operator ( |4~2l ), u* , with the bump of /o -field 
model, uo, and the bumps of the /i-field model, n^* and Ur, given on [A,-, Aq]. 
The connectivity u is as in ( |2.1| ) and / as in ( |2.10| ). The parameters are given 
as in the text, (b) The error sequence defined as in (|4.5|). 









u(x) 

h 

h+T 





















Figure 5 . The bump solution of the /-field constructed from u* , see Figj4ja). 
The parameters used are as in the text. 



For t G [0, r] we assume that there exist the excitation width A{t) such that a bump solution 
to /-field model, ua{x), satisfies 

UA{±A{t)) =t + h. 

Then ua (x) can be described by 




r{x,y)dydC 



using the representation ( |2.11[ ). 

Let Bbe a. real Banach space with partial ordering > defined by a cone E = {u ^ B\u{x) > 
0}. In this section we assume B = L2([0, r]). Then, if a bump of the /-field is given by ( |2.12| ) 
then A{t) € [At-, AqJ . The excitation width A satisfies the fixed point problem 

(5.2) A = AA, {AA){t) = A{t) + k{uA{A{t)) -t - h) , A: = const E M. 
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Theorem 5.1. The operator A is Frechet dijferentiable in L2[0, t] if f e W^'°°{R). 
Proof. Let us define the operator 

{GA){t) = UA{A{t)). 

We calculate the Frechet derivative of the operator GA. To do this we first compute its Gateaux 
derivative 

For any 5 G L2 [0, r] let us consider 

(5.3) G(A + X5) - G(A) = T piOhim, 

Jo 

where 

h{0= r{A{t) + X6{t),y)dy- r{A{t),y)dy. 



Making use of the Taylor expansion for r(A(t) + X6{t),y) as a function of A at A = we have 

A{€)+A5(0 / dr \ 



li(0= / ^r{Ait),y) + X—{A{t),y)6it) + o{X)jdy- r{A{t),y)dy 

r{A{t), y)dy + X6{t) / — (y, A(t))(iy + o(A). 



A(C) 



5x 



Plugging /i into ( |5.3| ) and making use of the mean value theorem we get the following formula 

(5.4) g{A, 6{t)) = j^^ p{OS{Or{A{t),A{0)d^ + 6{t) j^^ p{0^{A{t), A{0)d^. 

Hence, we arrive at the conclusion that the Gateaux derivative is a linear operator. In order to 
prove Frechet differentiability of the operator G we show, in accordance with [18], that g{-,S) : 
L2[0,t] — )• L2[0,t] is a continuous operator for all 6 G L2[0, r]. The proof of this fact is 
technical and we therefore formulate it as a separate lemma: 

Lemma 5.2. The operator g{-,d) : L2[0, t] — )• L2[0, r] is continuous for all 6 G L2[0, t]. 



Proof. We consider the first and the second integral of ( |5.4[ ) separately as the operators of A. 
Using the Cauchy-Schwarz and Minkowskii inequalities we show that these operators are con- 
tinuous and, thus, g{-, 6) is continuous as well, for any 6 G L2[0, r]. We present the proof for 
the first integral operator. The proof of continuity for the second term proceeds in the same way 
and is omitted. 
We introduce 

{FA){t)= r p{mOr{Ait),AiCm. 
Jo 

We obtain 

(FAi - FA2)(t) = r piOm (r(Ai(t), Ai(0) - r(Ai(t), A2(0)+ 
Jo 

+r(Ai(t), A2(0) - r{A2{t), A2(e))) d^ = h{t) + hit) 
where by the mean value theorem Ii and I2 can be defined as 

hit) = £ p{mO^{Al{t),A^mA,{0 - A2{0)d^ 

hit) = ^%(0<5(0|;(A2(t), A2(0)(Ai(i) - A2{t))dC 
with Afc = AfeAi + (1 — Afc)A2, for some A^ G [0,1], k = 1, 2. 
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We consider the norm of the difference. Using the Minkowskii inequality we get 



|FAi -FA2||L2[0,r] = 



< 



(hit) + l2{t)fdt 
1/2 



1/2 



+ 



IhiWdt 



< 

1/2 



Applying the Cauchy - Schwarz inequality to each of the terms we have 

\\FAi-FA2\\L,[0,r] < 



<[II 

\0 

(r r 
// 


'1,00 



p(05(0|j(Ai(t),Ai(0) 

p(05(0|^(Ai(t),Ai(0) 



2 r \ 

d^}\Ai{0-A2{0?d^dtj + 

2 r \ 
dC J\Ai{t)- A2it)\^d^dt] 
/ 



Since r G VF^'°°(M x M), p G -Loo(K), and G L2[0, r] the following estimate is valid 



p(0<5(0fj(Ai(t),Ai(0) 
p(05(0|^(A2(t),Ai(0) 



d^ < 



where 



, , , , 5r , 



dy 



Therefore, we get the following inequality 

\\FAi - FA2\\L,[o,r] < \C\V^\\^i - A2||l2[0 
from which the continuity of F follows. 



(RxR)}||'^lli2[0,r]- 



□ 



For convenience we make redefinition: ^(A, 6) = Obviously, the operator A is Frechet 
differentiable in any A G L2[0, r] and 



(5.5) 



A'^ = I + kG'^ 



□ 



We have the following lemma: 



Lemma 5.3. The operator A' ^5 > OforS > and A G [[At-,Ao] under Assumption^and^ 

and < k < 1/m, where 



(5.6) 



m = - min — (A(t), A(e)). 
t,ce[o,r] dx 



Proof. First of all, we notice that d^/dx G i?C(M) x BC{M.). Thus, there exists a finite mini- 
mum of I dx on the given set. Moreover, this minimum is negative according to Assumption 
[3] Therefore, m given by ( |5.6| ) is finite and positive, and the operator A'^ preserves positivity 
for < A; < □ 



Theorem 5.4. Let the conditions of Theorem 5.1 and Lemma 5.3 be satisfied. Then the operator 
A : [[At-,AoI| — 7- D C L2[0,t] has a fixed point in ^Ar, AoJ. Moreover, the sequences {A'^Ar} 
and {A'^Aq} converge to the smallest and greatest fixed point of the operator A, respectively. 
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Proof. The operator A is monotonically increasing. Indeed, we let A2 > Ai. Then AIS.2 — 
^Ai = ^'^(A2 - Ai) where A € [Ai, A2]] C [A^, Aq]]. Using Lemma [53] we conclude that 
A is monotone. 

The operator A is Frechet differentiable, and hence continuous in -L2[0, r] (see Lemma 5.1 1. 
Moreover, we have the following inequalities 

(^Ao)(f) = Ao + /c(nAo(Ao) - t - /i) = Ao + /c(uo(Ao) - t - /i) = Aq - /ct < Aq. 

and 

{A/^r){t) = A^ + A:(uA,(A^) - t - /i) = A^ + /c(u^(A^) - t - /i) = Aq + A;(t - t) > A^. 
Applying Theorem |3.1| tve complete the proof. □ 



Remark 5.5. We prove Theorem 5.4 for the case when D € ^2(0, r] but do not consider the 
case D € C[0, r]. T/ze co?ie of positive functions in C[0, r] nof regular Therefore additional 
assumptions on the operator A are required (see Theorem 3.1 ). We notice that A is not compact 
in C[0, r]. Indeed, the operator A is a Frechet differentiable with A'^ defined as in \5.5\ where 
A'^ is a sum of the identity operator and a compact operator, thus is not compact. Therefore, A is 
not a compact operator, see 11181 . The operator A does not seem to be condensing either, at least 
with respect to the Hausdorff measure. The case of more general measures of noncompactness 
11171 120I is not considered here. 



It remains to show that ua, where A is the fixed point of ( |5.2| ), is a bump. The definition of 
UA requires A{t) monotonically decreasing. We introduce the assumption. 

Assumption 3'. The partial derivative of^ with respect to x is negative for x = A{t), y = A(s) 
for t,s £ [0, r] and A{t) is a fixed point of ( |5.2[ ), i.e., 

— (A(t),A(s)) <0, Vt,sG [0,r]. 

Lemma 5.7. The fixed point A(i) of operator A is monotonically decreasing and differentiable 
on [0, r] under Assumption 3'. 



Thus, 



Proof. Since A(t) is a solution of the fixed point problem (5^1 then nA(A(t)) = t + h. We 
prove the lemma by direct differentiation of the last equality with respect to t. We obtain 

^%(0|^(A(t),A(0)A'(t)de = l. 

A'(t)= (jJ|%(0|^(A(t),A(e)K) '<0 

as — (A(t), A(^)) < by Assumption 3'. □ 

ox 

Assumption 3' requires an apriori knowledge of A{t) and therefore can not be checked before 
A(t) is found. Thus, we suggest to replace this assumption with the following one: 

Assumption 3". The partial derivative of^ with respect to x is negative for allx,y £ [Ar , Aq] , 
i.e., 

— (x,y) < 0, Vx,y G [A^,Ao]. 

The fulfillment of Assumption 3" implies that Assumption 3 and Assumption 3' are satisfied. 
In addition to Assumption 3' (or 3") we have the following requirement: 

Assumption 5'. The function ^ is such that 

(i) <h,\/x> A(0), y G [A(t), A(0)], 
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(ii) ^x,y)>h + T, Vxe[0,A(r)], y G [A(r), A(0)]. 



Theorem 5.10. Let A be a fixed point refereed to in Theorem \5.4\ Then U/\ defined (5.1 1 is 
a bump solution to under Assumptions 3" and 5'. 

Proof. We rewrite ( |5.1| ) as 

Jo 

Next, we make use of Assumption 5'. Keeping in mind the normalization property of p we show 
that 

ua{x) <h, Vx > A(0), y G [A(r), A(0)], 
UAix) >h + T, Vx e [0, A(t)], y G [A(r), A(0)]. 

□ 



Remark 5.11. For operator Tj we use Assumptions^^ and Assumptions^^ 3" and 5' for the 
operator A. Assumptions 3" and 5' are more restrictive than Assumptionsj^and^ Moreover 
Assumption^needs information about the fixed point A(t) which is a disadvantage. On the 
other hand, the operator Tj requires one extra assumption, Assumption^ 



5.1. Numerical example. Let /i, r and A,-, A* , Aq are chosen as in Section 4.1 Them, 



as we have mentioned before. Assumptions T]2 and [3] hold true. Hence, we can apply Theorem 



5.4 and obtain A(t). In Figj6];a) we illustrate the result of the iteration process. In Figj6jb) we 



have plotted the errors calculated as 
(5.7) e(n) =max|(^'^|Ao)(t) - (A"A^)(t)|, n = 1, 2, A^. 

X 

Similar as in (|4.5|), defines the identity operator, n corresponds to the iteration number, and 



N denotes the total number of iterations. In our calculations e{n) < 10 ^ for n > 13, and the 
minimal and maximal fixed points converges to each other. Thus, the fixed point is unique, see 



Corollary 3.2 We also observe that the fixed point A(t) belongs to [A* , AqJ , see Figl6ra) 




Figure 6. (a) A fixed point of pJ) , A(t), vertical hnes Aq, A,-, andA^* as 
they defined in Section 4.1 The connectivity function uj is given as in FigjTJa), 
/ is defined by ( |2.10 ) with p = 3, h = 0.1, r = 0.05. (b) The errors given as 
inKTl. 
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Knowing A{t) we liave ciiecked tliat Assumption 5' is fulfilled. Thus, by Theorem 5.10 we 
can obtain a bump solution to the /-field model ( |1.1[ ). 



We claim that that this bump coincides with the bump constructed in Section 4. 1 To demon- 
strate this, we found 6{t) that solves 

u*{S{t)) = t + h, te[0,T] 

with u* being the fixed point of the operator Tj, see Figj4ja). We calculate the relative error as 



(5.8) e = max 



m 

For our example we have obtained e = 2.5 x 10^^. We notice here that our implementation is 
not optimal and can be significantly improved. We do not pursue this problem here, however. 

6. Discussion 

We have introduced two iteration schemes for finding a bump solution in the /—field of 
the Wilson-Cowan model: The first scheme is based on the fixed point problem formulated by 
Kishimoto and Amari Q. The second one is described by the fixed point problem formulated 
for the interface dynamics of the bump. The latter formulation became possible due to the 
special representation of the firing rated function introduced by Coombes and Schmidt fSl. 

We have proved using the theory of monotone operators in Banach spaces that both iteration 
schemes converge under Assumption [T] and |2] From the iterative procedures we obtain the 
solution on the finite interval [A,-, Aq] (see Section|4]l, and on [A(0), A(t)] (see Sectionjs]). 
Then it has been shown that under some additional assumptions on the connectivity function uj 
this solution determines a bump of the /-field on M. 

The assumptions imposed for the first method (see Section |4]l differ from the ones imposed 
for the second method (see Section [5]l. The evident disadvantage of Assumption 3' and 5' is 
that they contain information about the output of the iteration procedure, A(t). Assumption 
3' can be substituted with the more restrictive Assumption 3", but not Assumption 5'. Thus, 
Assumption 5' can not be checked in advance. On the other hand, the set of assumptions for the 



fixed point problem ( |5.2| ) is in general less restrictive than the assumptions imposed on the fixed 
point method outlined in Section [4] All assumptions (except Assumption 5') are quite easy to 
check if uj{x) is given. 

We show by a numerical example that both iterative schemes converge to the same solution. 
Moreover, from numerics it follows that this solution is unique and stable. Indeed, the maximal 
and minimal fixed points turn out to be equal for any trials and choice of parameters. Thus, by 



Corollary 3.2 the fixed point is unique. Moreover, the constructed fixed point solution is stable 
since it is located between stable solutions of the /q- and /T--field, [7]. Notice that we have not 
given a mathematical verification of these observations. 

Notice also that we have looked for the bump solutions under the assumption A,- < Aq and 



( |2. 12| ). Thus, even if the constructed solution is unique, it does not necessarily mean that there 
are no other stable or unstable solution. However, the same type of reasoning as we performed 
here are no longer valid if we relax on these assumptions. Therefore we leave this problem for 
a future study. 
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